C ERN-TH/98-301 
hcp-ph/9809377 



Theoretical Status of B — > Xs'j Decays 



Matthias Neubert 
Theory Division, CERN, CH-1211 Geneva 23, Switzerland 



Abstract: 

We review the theoretical understanding of the branching ratio and 
photon-energy spectrum in B ^ Xsj decays at next-to-leading order 
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Standard Model, we obtain B{B Xs-f) = (3.29 ± 0.33) x IQ-^ for the 
total branching ratio, and B{B Xsj) = [2.85+0,^0) x 10"'' if a cut 
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We review the theoretical understanding of the branching ratio and photon-energy spectrum in B — > Xsf decays at next- to- 
leading order in QCD, including consistently the effects of Fermi motion. For the Standard Model, we obtain B(_B — > Xs'j) = 
(3.29 ± 0.33) X 10-4 for the total branching ratio, and B{B Xs-y) = (2.85l[5'^g) ^ ^ ^^-f" > 2-1 GeV is applied on 

the photon energy, as done in the recent CLEO analysis. A precise measurement of the photon spectrum would help reducing the 
theoretical uncertainty and yield important information on the momentum distribution of b quarks inside B mesons. 



1 Introduction 

About three years ago, the CLEO CoUaboration reported 
the first measurement of the inclusive biianching ratio for 
the radiative decays B — > Xg^, yielding^ B(i? Xg^) = 
(2.32 ±0.57 ±0.35) x llj-*. At this Conference, this value 
has been updated toH B(B ^ X^i) = (3.15 ± 0.35 ± 
0.32 ± 0.26) X lO"-*, where the first error is statistical, 
the second systematic, and the third accounts for model 
dependence. The ALEPH Collaboration has reported a 
measurement of the corresponding branching mtio for h 
hadrons produced at the Z resonance, yieldingB B(i?f, 
Xs-i) = (3.11±0.80±0.72)xl0-l Theoretically, the two 
numbers are expected to differ by at most a few percent. 
Taking the weighted average gives 

B(B ^ X^7) = (3.14 ±0.48) X 10-^ (1) 

Being rare processes mediated by loop diagrams, 
radiative decays of B mesons are potentially sensitive 
probes of New Physics beyond the Standard Model, pro- 
vided a reliable calculation of their branching ratio can 
be performed. The theoretical framework for such a 
calculation is set by the heavy-quark expansion, which 
predicts that to leading order in 1/mt, inclusive decay 
rates agree with the partQiij|nodel rates for the underly- 
ing decays of the h quarka l3. The leading nonperturba- 
tive corrections have been studied in detail and are well 
understood. The prediction for the B — > Xs^ branch- 
ing ratio suffers, however, from large perturbative un- 
certainties if only leading-order expressions for the Wil- 
son coefficients in the effective weak Hamiltonian are em- 
ployedu-t2l. Therefore, it was an important achievement 
when the full next-to-leading order calculation of the to- 
tal B Xsj branching ratio in the Standard Model was 
completed, conohiping consistently results_for the match- 
ing conditions Hi dlj, matrix elements Ej^tj, and anoma- 
lous dimensions E£l. The leading QED axicLelectroweak 
radiative corrections were included, too EZlEa. As a re- 
sult, the theoretical uncertainty was reduced to a level of 
about 10%, which is slightly less than the current exper- 
imental error. During the last year, the next-to-leading 



order analysis wa£L,extended to the cases, of two-Higgs- 
doublet modelsE3c3 and supersymmetryEJ, so that accu- 
rate theoretical predictions are now at hand also for the 
most popular extensions of the Standard Model. 

The fact that only the high-energy part of the pho- 
ton spectrum in i? — > Xsj decays is accessible exper- 
imentally introduces a significant additional theoretical 
uncertainty ta. For instance,, in the new CLEO analysis 
reported at this Conference cl a lower cut on the photon 
energy of 2.1 GeV in the laboratory is imposed, which 
eliminates about three quarters of phase space in this 
variable. After reviewing the theoretical status of calcu- 
lations of the total B Xg'y branching ratio, we thus 
discuss to what extent the effects of a photon-energy cut- 
off can be controlled theoretically, pointing out the im- 
portance of measurements of the photon spectrum for 
reducing the theoretical uncertainty in the extraction of 
the total branching ratio. 

2 B — > Xgj branching ratio 

The starting point in the analysis of B ^ Xs'y decays is 
the effective weak Hamiltonian 

4G 

Heft = y=- V;,Vtb C^{pib)OM . (2) 

The operators relevant to our discussion are 

lOTT^ 

O8 = f0s-La^.Gria6fl. (3) 

To an excellent approximation, the contributions of other 
operators can be neglected. The renormalization scale /if, 
in (||) is chosen of order ruh, so that all large logarithms 
reside in the Wilson coefficients Ci{fib)- For inclusive de- 
cays, the relevant hadronic matrix elements of the local 
operators Qi can be calculated using the heavy-quark ex- 
pansion dlD. The complete theoretical prediction for the 
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B decay rate at next-to-leading order in as bas 

been presented for the first time by Chetyrkin et al.ll3. 
It depends on a parameter S defined by tlie condition 
that the photon energy be above a threshold given by 
> (1 - The prediction for the B X^j 

branching ratio is usually obtained by normalizing the 
result for the corresponding decay rate to that for the 
semileptonic rate, thereby eliminating a strong depen- 
dence on the 6-quark mass. We define 
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where f{z) « 0.542 - 2.23(Vz - 0.29) is a phase-space 
factor depending on the quark-mass ratio z = (mc/rnhf' ■ 
The fine-structure constant a is renormaliaed at = 0, 
as is appropriate for real-photon emissionll3. The quan- 
tity -H'nlo(^) contains the next-to- leading order cor- 
rections. In terms of the theoretically calculable ra- 
tio _Rth('5), the B Xsj branching ratio is given by 
B{B ^ X,7) = 0.105A^sLi?th(^), where iVgL = B(B ^ 
Xc e i^)/10.5% is a normalization factor to be determined 
from experiment. To good approximation A^sl — 1- The 
current experimental situation of measurements of the 
semileptonic branching ratio of B mesons and their the- 
oretical interpretation are reviewed in Refs. 22,23. 

In the calculation of the quantity K^hoi^) we con- 
sistently work to first order in the small parameters as, 
^/rriQ and a/ as, the latter ratio being related to the 
leading- logarithmic QED corrections. The structure of 
the result is 

^NLoW^ J2 hjiS,tib)Re[CMC*ifXb)] , (5) 

i, 3 = 2,7,8 

where the Wilson coefficients Ci{nb) are expanded as 



^(0) «^ (1) . 

in as[fib) 



(6) 



The coefficients C^'''' (pb) are complicated functions of the 
ratio r] = asimw) / as{fJ.b) , which also depend on the val- 
ues Ci{mw) of the Wilson coefficients at the weak scale. 
In the Standard Model, these inital conditions arc func- 
tions of the mass ratio Xt — {mt/mwY- Whereas the 
leadingpOrder coefficients cf^ {/ib) are known since a long 
timel3"t3, the next-to- leading terms in (^), which must 
be kept for the coefficient C7 have been calculated 
only recently. The expression for ci^^\fj,b) can be found 
in eq. (21) of Ref. 16, and the result for c!f'^\fib) is given 
in eq. (11) of Ref. 18. 



Explicit expressions for the functions kij{S, ^b) in (|^) 
can be found, e.g., in Ref. 18, where we have corrected 
some mistakes in the formulae for real-gluon emission 
used by previous authors. (The corrected expressions are 
also given in the Erratum to Ref. 16.) Bound-state cor- 
rections enter the formulae for the coefficients fcy at ord£t 
l/rriQ and are proportional to the hadronic parameter cil 
A2 = \{m\, - m|) « 0.12 GeV^. Most of them char- 
acterize the spin-denendent interactions of the b quark 
inside the B mesouBTO. However, a peculiar feature of 
inclusive radiative decays is the appearance of a correc- 
tion proportional to l/m^ in the coefficient fey, which 
represents a long-dislaace contribution arising from (cc) 
intermediate states E3c3. 



2. 1 Definition of the total branching ratio 

The theoretical prediction for the B branching 
ratio diverges in the limit 5 —^ 1 because of a logarithmic 
soft-photon divergence of the b — > sgj subprocess, which 
would be canceled by an infrared divergence of the 0{a) 
corrections to the process b ^ sg. We have argued that 
a reasonable definition of the "total" branching ratio is 
to use an extrapolation to (5 = 1 starting from the region 
5 ~ 0.5-0.8, where the theoretical result-exhibits a weak, 
almost linear dependence on the cutoffta. The extrapo- 
lated value so defined agrees, to a good approximation, 
with the result obtained by taking S = 0.9, and hence 
we define the total branching ratio using this particular 
value of the cutoff. 

The theoretical result is sensitive to the values of 
various input parameters. For the (one- loop) quark pole 
masses we take mc/rrib = 0.29 ± 0.02, rub = (4.80 ± 
0.15) GeV, and nit — (175 ± 6) GeV. The correspond- 
ing uncertainties in the branching ratio are, respectively, 
j:5 0%, Tl-0%, and ±1.6%. We use the two-loop expres- 
sion for the running coupling as (fJ-) with the initial value 
as{mz) = 0.118 ± 0.003, which induces an uncertainty 
of ±2.7%. For the ratio of the CKM parameters in (|) 
we take the value jt^* Vtfc|/| Vcfc| = 0.976 ± 0.010 obtained 
from a global analysis of the unitarity triangle E3. This 
gives an uncertainty of ±2.1%. Finally, we include an 
uncertainty of ±2.0% to account for next-to-leading elec- 
troweak radiative corrections til. The theoretical uncer- 
tainty arising from the variation of the renormalization 
scale will be addressed below. We find an uncertainty of 
±6.3% from the variation of the scale /if,, and of !ti'5% 
from the variation of the scale /If, entering the expres- 
sion for the semileptonic decay rate in the denominator 
in (^. Adding the different errors in quadrature gives 
a total uncertainty of l;^g% (adding them linearly would 
lead to the more conservative estimate of ^22'^)- Fo^' the 
total B — > Xsj branching ratio in the Standard Model 
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we obtain 

B{B -> Xsj) = (3.29 ± 0.33) x IO^'^TVsl , (7) 
in good agreement with the experimental value in (Q). 

2.2 Sensitivity to New Physics 

Possible New Physics contributions would enter the the- 
oretical prediction for the B Xsj branching ratio 
through non-standard values of the Wilson cocfhcicnts of 
the dipole operators O7 and Og at the weak scale rnvv- To 
explore the sensitivity to such effects, we normalize these 
coefficients to their values in the Standard Model and in- 
troduce the ratios = Ci{mw)/Cf^{mw) with i = 7,8. 
In the presence of New Physics, the parameters ^7 and 
may take (even complex) values different from 1. Simi- 
larly, New Physics may induce dipole operators with op- 
posite chirality to that of the Standard Model, i.e. opera- 
tors with right-handed light-quark fields. If we denote by 

and the Wilson coefficients of these new opera- 
tors, expression (||) can be modified to include their con- 
tributions by simply replacing dC* — > CiC* + C^Cf* 
everywhere. We thus define two additional parameters 

= Cf-{mw)/Cf^{mw) with i = 7,8, which vanish 
in the Standard Model. Since the dipole operators only 
contribute to rare flavour-changing neutral current pro- 
cesses, there are at present rather weak constraints on 
the values of these parameters!^. On the other hand, 
we assume that the coefficient C2 of the current-current 
operator O2 takes its standard value, and that there is 
no similar operator containing right-handed quark fields. 
With these definitions, the B — > Xg^ branching ratio can 
be decomposed aslla 
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B(B^X,7)| 



iS-,>(l-5)iS»-- 
2 I lt-R|2\ 



A^SL 

= B22{5) +i?77(5)(|C7p + +B,,{5){U^ + |4«f ) 

+ B27(5)Re(^7) + ^28('5)Re(C8) 

+ B,^{5) [Re(e7e8) + Re(C7^ef )] . (8) 

In Table 0, the values of the components Bij{S) ob- 
tained with 5 = 0.9 are shown for different choices of 
the renormalization scale. Assuming TVgL = 1, the Stan- 
dard Model result for the total branching ratio is given by 
BsM — X^ij ^ij ■ The most important contributions are 
-827 and -B22 followed by ^77. The smallness of the re- 
maining terms shows that there is little sensitivity to the 
coefficient Csimw) of the chromo-magnetic dipole oper- 
ator. Once the parameters and are calculated in a 
given New Physics scenario, the result for the B Xgj 
branching ratio can be derived using the numbers shown 
in the table. For the remainder of this talk, however, we 
assume the validity of the Standard Model. 



Table 1: Values of the coefficients Bij (for 5 = 0.9) in units of 
10~*, for different choices of fii,. The coefficient Sgg = 0.015 in all 
three cases 



mb/2 
nib 
2mb 



B27 



B- 



22 



B 



77 



B 



28 



B 



78 



1.265 1.321 
1.395 1.258 
1.517 1.023 



0.335 
0.382 
0.428 



0.179 
0.161 
0.132 



0.074 3.188 
0.083 3.293 
0.092 3.206 



2.3 Perturhative uncertainties 

The components Bij{5) in (^ are formally independent 
of the renormalization scale /i^. Their residual scale de- 
pendence results only from the truncation of perturba- 
tion theory at next-to-leading order and is conventionally 
taken as an estimate of higher-order corrections. Typi- 
cally, the different components vary by amounts of order 
10-20% as ^b varies between m;,/2 and 2mb. The good 
stability is a result of the cancelation of the scale depen- 
dence between Wilson coefficients and matrix elements 
achieved by a next-toJcading-order calculation. 

Previous authorsEjlij'E3E3 who estimated the /ib de- 
pendence of the total B Xs^ branching ratio in the 
Standard Model found a more striking improvement over 
the leading-order result, namely a variation of only ^^^'9% 



as compared with ^20^ leading order. 



3.2' 

20 /u leaning uiu.ci. Howcvcr, the 
apparent excellent stability observed at next-to-leading 
order is largely due to an accidental cancelation between 
different contributions to the branching ratio. A look 
at Table |^ shows that the residual scale dependence of 
the individual contributions Bij is much larger than that 
of their sum, which determines the total branching ra- 
tio in the Standard Model. Note, in particular, the al- 
most perfect cancelation of the scale dependence between 
-B27 and i?22, which is accidental since the magnitude 
of i?27 depends on the top-quark mass, whereas -B22 is 
independent of mt. In such a situation, the apparent 
weak scale dependence of the sum of all contributions is 
not a good measure of higher-order corrections. Indeed, 
higher-order corrections must stabilize the different com- 
ponents Bij individually, not only their sum. The varia- 
tion of the individual components as a function of /i;, thus 
provides a more conservative estimate of the truncation 
error than does the variation of the total branching ratio. 
For each component, we estimate the truncation error by 
taking one half of the maximum variation obtained by 
varying /if, between mb/2 and 2'mb. The truncation er- 
ror of the sum is then obtained by adding the individual 
errors in quadrature. In this way, we find a total trun- 
cation error of ±6.3%, which is more than a factor of 2 
larger than the estimates obtained by previous authors. 
An even larger truncation error could be justified given 
that the choice of the range of variation of fib is ad hoc. 
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and that the scale dependence of the various components 
is not symmetric around the point fib = rnt. 



3 Partially integrated branching ratio and 
photon-energy spectrum 

Whereas the explicit power corrections to the functions 
kij are small, an important nonperturbative effect not 
included so far is the motion of the b quark inside the 
B meson caused by its soft interactions with the light 
constituents. It leads to a modification of the photon- 
energy spectrum, which muatJae taken into account if 
a realistic cutoff is imposed E2c3. This so-called "Fermi 
motion" can be included in the heavy-quark expansion 
by resumming an infinite set of leading-twist contribu- 
tions into a shape function _F(fc+), which governs the 
light-cone momenluJiiLjlistribution of the heavy quark in- 
side the B mesonEjT3. The physical decay distributions 
are obtained from a convolution of parton-model spectra 
with this function. In the process, phase-space bound- 
aries defined by parton kinematics are transformed into 
the proper physical boundaries defined by hadron kine- 
matics. 

The shape function is a universal characteristic of 
the B meson governing the inclusive decay spectra in 
processes with massless partons in the final state, such 
as i? — > and B i v. However, this function 

does not describe in an accurate way .the .distributions 
in decays into massive partons such ascJo B Xciv. 
Unfortunately, therefore, the shape function cannot be 
determined using the lepton spectrum in semileptonic de- 
cays, for which high-precision data exist. On the other 
hand, there is some useful theoretical information on the 
moments An = (fc") of the function F{k+), which are m- 
lated to the forward matrix elements of local operatorstJ. 
In particular, Ai = vanishes by the equations of mo- 
tion (this condition defines the heavy-quark mass), and 
A2 — \\J?^ is related to the kinetic energy of the heavy 
quark inside the B meson. For our purposes it is sufficient 
to adopt the simple form F{k+) = N (1 — x)°e(^+"''^ with 
X — /c+/A < 1, where A = — m^. This ansatz is such 
that Ai = by construction. The parameter a can be 
related to the second moment, yielding fi^ = 3A^/(l-|-a). 
Thus, the &-quark mass (or A) and the quantity (or 
a) are the two parameters of our function. We take 
rub = 4.8 GeV and fi^. = 0.3 GeV^ as reference values, 
in which case a ~ 1.29. 

Let us denote by Bfj{Sp) the various components in 
(||) calculated in the parton model, where the cutoff is 
defined by the condition that > ^(1 — <5p)mf,. Then 
the corresponding physical quantities Bij{5) with 5 de- 
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Figure 1: Partially integrated B — ► X37 branching ratio (upper 
plots) and photon-energy spectrum (lower plots) for various choices 
of shape-function parameters 

fined such that E-^ > ^{1 — S)mB are given bylllS 

mB—mi, 
m_B(l — (5) — Tjii, 

(9) 

This relation is such that -Bij(l) — Bfj{l), implying that 
the total branching ratio is not affected by Fermi motion. 
The effect is, however, important for realistic values of 
the cutoff S. 

As an illustration of the sensitivity of our results to 
the parameters of the shape function, the upper plots 
in Figure |l| show the predictions for the partially inte- 
grated B — > Xsj branching ratio as a function of the 
energy cutoff E'""" ~ '^(^ ~ 5)mB- In the left-hand plot 
we vary keeping the ratio /iJ/A^ fixed. The gray line 
shows the result obtained using the same parameters as 
for the solid line but a different functional form given by 
F(fc+) = A^(l — xYe'^'^'^^^^ . For comparison, we show 
the data point B(B X,'^) = (2.04zL0.47) x 10"'' ob- 
tained in the original CLEO analysis EJ with a cutoff at 
2.2 GeV. In the right-hand plot, we keep rrib = 4.8 GeV 
fixed and compare the parton-model result (gray dot- 
ted curve) with the results corrected for Fermi motion 
using different values for the parameter /ij. This plot 
illustrates how Fermi motion fills the gap between the 
parton-model endpoint at mi,/ 2 and the physical end- 
point at mB/2. (The true endpoint is actually located 
at [rrig — {ttik + ■mT^Y]/2mB ~ 2.60 GeV, i.e. slightly 
below TUB 12 K, 2.64 GeV.) Comparing the two plots, it 
is evident that the uncertainty due to the value of the 
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Table 2: Partially integrated B — » Xs7 branching ratio for different 
values of the cutoff on the photon energy 



7, mill 


B{B ^ X,7) [10-*A^sl] 


2.2 GeV 
2.1 GeV 
2.0 GeV 


2.56 ±0.26^0,11 
2.85 ± 0.29 

3.01 ± 0.30 tHs 



6-quark mass is the dominant one. Variations of the 
parameter /i^ have a much smaher effect on the par- 
tiaUy integrated branching ratio, and also the sensitiv- 
ity to the functional form adopted for the shape function 
turns out to be small. This behaviour is a consequence 
of global quark-hadron duality, which ensures that even 
partially integrated quantities are rather insensitive to 
bound-state effects. The strong remaining dependence 
on the fe-quark mass is simply due to the transformation 
by Fermi motion of phase-space boundaries from parton 
to hadron kinematics. 

Taking the three curves in the left-hand upper plot 
for a representative range of parameters and applying a 
small correction for the boost from the B rest frame to 
the laboratory as appropriate for the CLEO analysis, we 
show in Table |2| the predictions for the partially inte- 
grated branching ratio for three different values of the 
cutoff on the photon energy. The first error refers to the 
dependence on the various input parameters discussed 
previously, and the second one accounts for the uncer- 
tainty associated with the description of Fermi motion. 
In general, this second error can be reduced in two ways. 
The first possibility is to lower the cutoff on the photon 
energy. A first step in this direction has already been 
takeiL in the new CLEO analysis reported at this Confer- 
ence □, in which the cutoff has been lowered from 2.2 to 
2.1 GeV. If a value as low as 2 GeV could be achieved, the 
theoretical predictions would become rather insensitive 
to the parameters of the shape function. To what extent 
this will be possible in future experiments will depend 
on their capability to reject the background of photons 
from other decays. The second possibility is that future 
high-precision measurements of the photon spectrum will 
make it possible to adjust the parameters of the shape 
function from a fit to the data. For the purpose of illus- 
tration, the photon spectra corresponding to the various 
parameter sets are shown in the lower plots in Figure ^ 
Such a determination of the shape-function parameters 
from B —^ Xg"/ decays would not only help to reduce the 
theoretical uncertainty in the determination of the total 
branching ratio, but would also enable us to predict the 
lepton. spectrum in _B — > £v in a, model-independent 
wayL3. This may help to reduce the theoretical uncer- 
tainty in the value of \Vub\- A detailed analysis of the 



photon spectrum will therefore be an important aspect 
in future analyses of inclusive radiative B decays. 

4 Conclusions 

The inclusive radiative decays B X^^ play a key role 
in testing the Standard Model and probing the structure 
of possible New Physics. A reliable theoretical calcula- 
tion of their branching ratio can be performed using the 
operator product expansion for inclusive decays of heavy 
hadrons combined with the twist expansion for the de- 
scription of decay distributions near phase-space bound- 
aries. To leading order in 1/mf, the decay rate agrees 
with the parton-model rate for the underlying quark de- 
cay b Xsj- With the completion of the next-to-leading 
order calculation of the Wilson coefficients and matrix 
elements of the operators in the effective weak Hamilto- 
nian, the perturbative uncertainties in the calculation of 
this process have been reduced to a level of about 10%. 
Bound-state corrections to the total decay rate are sup- 
pressed by powers of 1/mf, and can be controlled in a 
systematic way. 

A more important effect is the Fermi motion of the 
heavy quark inside the meson, which is responsible for 
the characteristic shape of the photon-energy spectrum 
in _B — > Xgj decays. It leads to the main theoretical un- 
certainty in the calculation of the branching ratio if a re- 
striction to the high-energy part of the photon spectrum 
is imposed. Fermi motion is naturally incorporated in 
the heavy-quark expansion by resumming an infinite set 
of leading-twist operators into a non-perturbative shape 
function. The main theoretical uncertainty in this de- 
scription lies in the value of the 6-quark mass. Other 
features associated with the detailed functional form of 
the shape function play a minor role. The value of nib 
and other shape-function parameters could, in principle, 
be extracted from a precise measurement of the photon- 
energy spectrum, but also gross features of this spectrum 
such as the amrage photon energy would provide valuable 
information E3. For completeness, we note that besides 
the photon-energy spectrum also the invariant hadronic 
mass distribution in radiative B decays can be studied. 
Investigating the pattern of individual hadron resonances 
contributing to the spectrum, one can motivate a simple 
description of the hadronic mass spectrum with only a 
single parameter, the B — > branching ratio, to be 
determined from experiment £3. 

Possible New Physics contributions would enter the 
theoretical predictions for the B Xs^ branching ra- 
tio and photon spectrum through the values of param- 
eters and ^j^, which are defined in terms of values of 
Wilson coefficients at the scale raw- This formalism al- 
lows one to account for non-standard contributions to the 
magnetic and chromo-magnetic dipole operators, as well 
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as operators with right-handed light-quark fields. Quite 



generally, New 
photon spectrum 



tfsics would not affect the shape of the 
3 but could change the total branching 
ratio by a considerable amount. This implies that the 
analysis of the photon-energy and hadronic mass spec- 
tra, which is crucial for the experimental determination 
of the total branching ratio, can be performed without 
assuming the correctness of the Standard Model. On the 
other hand, the total branching ratio will provide a pow- 
erful constraint on the structure of New Physics beyond 
the Standard Model as experimental data become more 
precise. 
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